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ABSTRACT 

We explore theoretical models of the ionization ratios of the Li-like absorbers 
N V, O VI, and C IV, in the Galactic halo. These ions are believed to form in 
nonequilibrium processes such as shocks, evaporative interfaces, or rapidly cooling 
gas, all of which trace the dynamics of the interstellar medium. As a useful new 
diagnostic, we focus on velocity-resolved signatures of several common physical 
structures: (1) a cooling Galactic fountain flow that rises, cools, and recombines 
as it returns to the disk; (2) shocks moving toward the observer; (3) a conductive 
interface with the observer located in the hotter gas. This last geometry occurs 
with the solar system inside a hot bubble, or when one looks out through the 
fragmenting top shell of our local bubble blown into the halo as part of the 
Galactic fountain. In Paper II, these models are compared to ionization-ratio 
data from FUSE and Hubble Space Telescope. 

Subject headings: Galaxy: halo — ISM: structure — ultraviolet: ISM 



1. Introduction 

The nature and dynamics of the inter- 
stellar medium (ISM) of galaxies deter- 
mines how the energy and matter released 
by stars are redistributed through the uni- 
verse. It is thus critical to understand the 
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ISM, and in particular that of the Milky 
Way, which can be observed with greater 
sensitivity and resolution than other galax- 
ies. The ISM in the disk of our galaxy 
consists of several phases, from dense and 
cold molecular gas to hot and fully ionized. 
Stars heat and disperse the dense clouds 
in which they form, and that hot gas cools 
and recombines eventually to complete the 
cycle and form more stars. It is particu- 
larly interesting to consider the interface 
between this multiphase medium and in- 
tergalactic space, where hot gas is released 
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from supernovae to several kiloparsecs in 
altitude, forming a hot diffuse medium first 
proposed by Spitzer (1956) as the Galactic 
corona. 

The Galactic corona or halo is almost 
certainly a dynamic object. It is difficult to 
construct static halo models, because mod- 
els supported by thermal pressure are ther- 
mally unstable. If conduction is sufficiently 
important to stabilize small-scale instabil- 
ities, then the entire (nearly isothermal) 
halo is unstable to collapse or expulsion 
ClS db wind (see Bregman 1980; Field 1965, 
for stability arguments). Cosmic-ray sup- 
ported static halos have been proposed 
(e.g., Boulares & Cox 1990), but there 
are considerable uncertainties as to how 
cosmic rays are confined by the Galaxy 
(the Galactic magnetic field topology in 
particular), and these models do not ex- 
plain the high and intermediate velocity 
clouds. Considering these things, Shapiro 
& Field (1976) first proposed a "Galac- 
tic fountain" of supernova-heated gas ris- 
ing buoyantly above the disk until it cools 
and falls back to the disk. Smooth Galac- 
tic fountain models have been constructed 
by several authors. In particular, Breg- 
man (1980) discussed the issues of radial 
flow in the Galactic gravitational field for 
a supersonic hot flow which reaches sev- 
eral kiloparsecs height. More recent mod- 
els (Houck & Bregman 1990; Breitschw- 
erdt et al. 1993) consist of transsonic flows, 
which require cooler initial temperatures 
and only rise to 1-2 kpc. These models 
can help to explain the population of in- 
termediate velocity halo clouds, believed 
to exist at those lower heights and veloc- 
ities compared to the more distant high 
velocity clouds (e.g., Wakker et al. 1998). 



The cooling layer certainly fragments into 
small clouds by Rayleigh- Taylor instability 
even if no other inhomogeneities exist (e.g.. 
Berry et al. 1998). 

Superbubbles, worms, and shell-like 
structures have been observed in (e.g., 
Heilcs 1984) the upper disk/low halo (few 
hundred parsecs altitude), but it is diffi- 
cult to ascertain whether the hot gas is 
being expelled into the halo. A related is- 
sue is the flUing factor of hot gas in the 
disk, which determines the rate at which 
hot gas can escape into the halo. Esti- 
mates for that filling factor and for the 
rate of reheating by halo supernovae are 
based on the evolution of supernova rem- 
nants in smooth media, but our current 
understanding of the ISM is increasingly 
inhomogeneous and dynamic, so those ar- 
guments are probably of limited value 
(Kahn 1998). Recently, the combination 
of good observations, sophisticated numer- 
ical models, and sufficient comprehension 
of the Galactic fountain may be begin- 
ning to allow identification of fountain-like 
rising structures (de Avillez & Mac Low 
2001). 

The dynamics of interstellar gas in gen- 
eral and specifically in the Galactic halo 
may be best understood by studying the 
hot phase of the ISM (coronal gas at mil- 
lions of degrees), and gas at temperatures 
~10^ K intermediate to the hot phase and 
cooler phases (10^ K). Of the different in- 
terstellar gas phases, the coronal gas is 
most directly linked to the main sources 
of energy in the ISM, supernovae and stel- 
lar winds. Slightly cooler gas is most 
closely linked to transient and dynamical 
processes. This gas is typically short-lived 
because the cooling time is short at at 
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10^ K. The lithium-like ions of common 

metals, in particular O VI, N V, and C IV, 
are sensitive tracers of interstellar gas at 
several times 10^ K. The resonance absorp- 
tion lines of these ions are observable with 
ultraviolet spectrographs on the Far Ultra- 
violet Spectroscopic Explorer (FUSE), and 
the Space Telescope Imaging Spectrograph 
(STIS) and Goddard High Resolution Spec- 
trograph (GHRS) aboard the Hubble Space 
Telescope (HST). 

Section 2 describes previous models of 
high-ion column densities. Section 3 de- 
scribes new models of the dynamical sig- 
natures of the Li-like ions and their inter- 
pretation. We summarize in § 4. 

2. Previous models of Li-like ions in 
the dynamic ISM 

The ions O VI, N V, and C IV are pre- 
dicted in models of transient phenomena 
such as shocks or interfaces between dif- 
ferent temperature gas, with conductive 
heating, turbulent mixing, or rapid radia- 
tive cooling (see Spitzer 1996, for a sum- 
mary). Notable, especially for O VI, is 
the importance of collisional ionization; 
by contrast. Si IV is produced in such 
transient phonomena, but is also com- 
monly produced by photoionization in the 
Galaxy. The column densities and col- 
umn density ratios of these high ions can 
be used to determine which physical sce- 
nario is predominant in the ISM. The 
distributions of C IV and N V in the 
Galaxy are (0.6-3) x 10^'^ cm"^ kpc"^ and 
(0.5-1) X 10^'^ cm"^ kpc"^ respectively (Sav- 
age & Massa 1987). The distribution of 
O VI is 1-2 regions per kpc, each with a 
column density of (2-5) x 10^^ cm"^. These 
statistics from Shelton & Cox (1994) came 



from reanalysis of the Copernicus dataset 
presented by Jenkins (1978), who found a 
higher frequency of lower column-density 
features. The range in the column den- 
sity ratio log[N(C IV)/N(0 VI)] has been 
previously measured at -1.5 to -0.5 in the 
disk and -0.5 to +0.5 in the halo (Spitzer 
1996). The values of log[N(N V)/N(0 VI)] 
reported by Sembach et al. (2001b) are - 
0.4 to -0.9 in the halo, but these represent 
a small heterogeneous sample of the dif- 
ferent lines of sight analyzed by different 
authors. 

The simplest model of gas containing 
Li-like ions is moderately hot gas in colli- 
sional ionization equilibrium (CIE); C IV, 
N V, and O VI peak in ionization fraction 
at temperatures of ~1, 2, and 3x10^ K, 
respectively (e.g., Shull & Van Steenberg 
1982; Sutherland & Dopita 1993). The ra- 
tio between any two species implies a CIE 
temperature. If derived from two differ- 
ent ion ratios, this CIE temperature does 
not agree for most real observations of the 
Galactic halo, indicating that the halo is 
probably not in collisional ionization equi- 
librium. This is seen in Figure 1, in which 
the CIE line ratios lie far to the lower left 
of more realistic, nonequilibrium models. 
Although there is some disagreement as 
to how well atomic abundances are known 
(Holweger 2001), recent changes in the so- 
lar abundances of carbon and oxygen are 
insufficient to bring these two derived tem- 
peratures into agreement, and in fact the 
latest C/0 (AUende Prieto et al. 2002) is 
very close to older values. 

It is important to understand what the 
predicted ion column density ratios arc in 
more dynamic physical models, beginning 
with conductive interfaces. Ballet et al. 
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Fig. 1. — Line ratios for various models of high- 
ion production. Tlie data plotted are from mod- 
els in the literature: radiative cooling of galac- 
tic fountain gas (triangles, Shapiro & Benjamin 
1993; Benjamin & Shapiro 1993), turbulent mixing 
layers (stars, Slavin et al. 1993), and conductive 
heating and evaporation of spherical, Bohringer 
& Hartquist (1987) and planar, Borkowski et al. 
(1990) clouds ("x"s), and in cooling supernova 
remnant shells (pluses, Slavin & Cox 1993; Shelton 
1998). Also shown are the ion ratios for hot gas 
in coUisional ionization equilibrium (Sutherland & 
Dopita 1993, CIE, solid line,). 



(1986) modeled the evaporation of evap- 
orating spherical clouds with nonequilib- 
rium ionization. They presented results for 
a small (i? = 5 pc) isolated cloud in a Mc- 
Kee & Ostriker (1977) type model of cool 
cloudlets in a hot medium of large filling 
factor. Relative to calculations of evapo- 
ration with collisional ionization equilib- 
rium, the total column densities of the 
Li-like ions is increased, as is the dis- 
tance from the cloud at which those species 
reach their maximum abundance. Both 
effects result from delayed ionization (up 
to the Li-like ionization stage but espe- 
cially up to the He-like stage) as the gas 
is conductively heated. Column densi- 
ties for all modeled structures are listed 
in Table 1, and column density ratios are 
given in Table 2. Bohringer & Hartquist 

(1987) also modeled evaporating spherical 
clouds with nonequilibrium ionization, but 
self- consistently included the effects of ra- 
diative cooling on the temperature struc- 
ture. They found somewhat lower ioniza- 
tion states than Ballet et al. (1986), with 
higher N V/0 VI and C IV/0 VI. 

Borkowski et al. (1990) presented more 
sophisticated models of a plane-parallel in- 
terface between hot and cold gas. They 
modeled the time-dependent evolution 
through an evaporative, steady-state, and 
recondensation phase. Numbers presented 
in Tables 1 and 2 are for the longest, 
steady-state phase, but much higher values 
of N(N V)/N(0 VI) and N(C IV)/N(0 VI) 
are obtained during the initial evapora- 
tion, as the species with lower ionization 
potentials ionize up to the Li-like stages 
faster. These authors also included a mag- 
netic field, which can suppress thermal 
conduction if oriented parallel to the in- 
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Table 1 
ll-like ion column densities. 



Physical Situation 



Si IV 
(1012 



CIV 



cm 



') (10 



12 



cm 



N V 
2) (10 



12 cm-2^ 



O VI 
(10^2 



cm • 



Reference 



observed sightlincs, kpc'^ 


2-10'^ 


6-30^^ 


5 10=* 


10 50^ 




evaporating cloudlet 




1.2-1.5 


0.5-0.6 


9-12 


Ballet et al. (1986) 


evaporating cloudlet 


.10-.14 


2.7-3.8 


1.0-1.2 


12-14 


Bohringer & Hartquist (1987) 


planar conduction front 


.10 .16 


1.6 3.2 


.6 1.0 


8 10 


Borkowski et al. (1990) 


planar conduction front 


.029-.097 


.89-2.7 


.40-1.0 


6.7-14 


Slavin (1989) 


stellar wind bubble 


.21-.25 


3.3-4.0 


1.3-1.6 


21-25 


Weaver et al. (1977) 


SNR bubble 


.4-.6 


6.3-10 


3.2 5.0 


40-79 


Slavin & Cox (1992) 


SNR bubble 


-0.52 


-7.8 


-3.6 


-47 


Slavin & Cox (1993) 


halo SXR 1hi1)1)1o 




8 15 


3.4 7.9 


3o I-IO 


Shclton (1998) 


4 Mq cooling 


3.3 6.4 


43 79 


28 36 


580-600 


Edgar & Chevalier (1986) 


40 pc cooling cloud 


-25 


-50 


-13 


-200 


Benjamin & Shapiro (1993) 


turbulent mixing layer 


.0010-.47 


.025-6.8 


.0022-.32 


.017-.81 


Slavin et al. (1993) 


white dwarfs 


1.4-4.4 


25-77 


3.7-12 


5.6-20 


Dupree & Raymond (1983) 



Note. — Predicted column densities vary with physical conditions, for example the temperatures of the 
hot and cold media in interface models, and an indicative range of column densities is given for each physical 
model. 

^Savage & Massa (1987). 

I'Shelton & Cox (1994); Jenkins (1978). 
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Table 2 

Ll-LIKE ION RATIOS. 



Physical Situation 



log 



N(C IV) 
N(0 VI) 



log 



N(N V) 
N(0 VI) 



Reference 



observed halo sightlines 


-0.5 to +0.5 


-0.4 t -0.9 


Spitzer (1996); Scmbach et al. (2001b) 


evaporating cloudlet 


-0.95 


-1.35 


Ballet et al. (1986) 


evaporating cloudlet 


-0.60 


-1.07 


Bohringer & Hartquist (1987) 


planar conduction front 


-0.4 to -0.8 


-0.9 to -1.1 


Borkowski et al. (1990) 


planar conduction front 


-0.7 to -0.9 


-1.1 to -1.2 


Slavin (1989) 


stellar wind bubble 


-0.8 


-1.2 


Weaver ct al. (1977) 


SNR bubble 


-0.8 to -0.9 


-1.1 to -1.2 


Slavin & Cox (1992) 


SNR bubble 


~-0.8 


~-l.l 


Slavin & Cox (1993) 


halo SNR bubble 


-0.6 to -1.1 


-0.9 to -1.2 


Shchon (1998) 


4 Mq cooling 


-0.88 to -1.1 


-1.2 to -1.3 


Edgar & Chcvahcr (1986) 


40 pc cooling cloud 


-0.4 to -0.5^^ 


-1.2 to -1.4 


Benjamin & Shapiro (1993) 


turbulent mixing layers 


-1-0.1 to +0.9 


-0.4 to -0.9 


Slavin et al. (1993) 


WD Stromgren spheres 


+0.6 to +0.7 


-0.18 to -0.22 


Dupree & Raymond (1983) 



'^Much lower C IV/0 VI ratios (~ —1.7 in log) are present in the 10^ K gas before it rapidly cools, but 
the cooler, overionized stage is much longer lived. 



terface. The total column densities of 
all species are lower in this case, but the 
N(C IV)/N(0 VI) and N(N V)/N(0 VI) 
ratios also decrease, because the evapo- 
rated gas is hotter, and the overall ioniza- 
tion state of the front is increased. This ef- 
fect was confirmed quantitatively by Slavin 
(1989). 

Weaver et al. (1977) modeled a conduc- 
tive interface of a different kind, at the 
shell of a wind-blown bubble around a mas- 
sive star. The hot gas on the bubble inte- 
rior evaporates the cooler shell that forms 
around the structure after the swept-up 
gas is sufficiently dense to cool. Slavin 
& Cox (1992) performed a similar calcu- 
lation for a SNR bubble and found higher 
total column densities, as expected for that 
more energetic phenomenon, but nearly 
the same ion ratios, as expected for the 



same physics. Slavin & Cox (1993) ex- 
plored the same type of models in more de- 
tail and obtained similar results, with scal- 
ing laws. Finally, Shelton (1998) extended 
these SNR models and applied them to a 
halo with a lower ratio of thermal to non- 
thermal pressure. She found a variation of 
the column densities and their ratios with 
time as the remnant evolves, with a signif- 
icant decline in the ionization state after 
10^ years, as the shell slows down and the 
interior cools. 

Edgar & Chevalier (1986) modeled the 
nonequilibrium ionization of cooling gas 
and scaled their results to a Galactic foun- 
tain with 4 Mq yr"^ of cooling matter. 
More detailed models of cooling gas in- 
corporating radiative self-ionization (Ben- 
jamin & Shapiro 1993; Shapiro & Ben- 
jamin 1991; Shapiro & Benjamin 1993) 
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have a similar range of Li-like ion produc- 
tion. Slavin et al. (1993) modeled turbu- 
lent mixing layers between million degree 
gas and cooler (10^-10^ K) gas. O VI, N V, 
and C IV are produced both in the coohng, 
overionized hot gas, and in the heating, 
underionized cool gas. These models are 
physically interesting because the topology 
of the interstellar medium includes many 
interfaces between hot and cold gas mov- 
ing at substantial relative velocities. The 
models can produce a large range of ion 
column density ratios and generally pro- 
duce smaller total column densities per in- 
terface than the other models. However, a 
number of the model parameters are un- 
constrained by the current state of the ob- 
servations. 

Although C IV and Si IV can be pro- 
duced by photoionization near hot stars 
and somewhat in the ISM by ambient stel- 
lar ionizing radiation and the soft X-ray 
background (e.g., Cowie et al. 1981; Black 
et al. 1980), O VI and N V are almost def- 
initely coUisionally ionized in the Galac- 
tic disk and halo. Nevertheless, there 
are some models that produce these ions 
by photoionization near a sufficiently hot 
source. Dupree & Raymond (1983) found 
that hydrogen-rich white dwarfs can pro- 
duce some O VI in their Stromgren spheres 
and in diffuse regions if the circumstellar 
medium is patchy. The Li-like ion column 
densities depend on impact parameter of 
a line of sight with the Stromgren sphere, 
and also increase rapidly with increasing 
average ISM density (values in Table 1 are 
for Uh = 0.1-1.0 cm"^). 

Line ratios for model production mech- 
anisms for Li-like ions are plotted in Figure 
1. Clearly, measurements of these line ra- 



tios can help to distinguish between differ- 
ent models. Additional detail may even be 
possible - for example, in the conductive 
interface models of Bohringer & Hartquist 
(1987) and Borkowski et al. (1990) ("x"s 
in Figure 1), models with increasing mag- 
netic field suppressing the conduction are 
located farther to the lower left of the dia- 
gram. 

It should be noted that the interstel- 
lar abundances of the atoms whose ions 
are being studied here are not well de- 
termined. The models presented above 
are for solar abundances, reported in Cox 
(2000), mostly from the work of Grevesse 
et al. (1996), measured in meteorites and 
in the solar atmosphere. This technique 
is particularly difficult for C, N, and O. 
The models given later in this paper build 
on the models above, and so the implicit 
metaUicities are those of Grevesse et al. 
(1996). There is the additional uncertainty 
of whether solar abundances are appropri- 
ate for the general interstellar medium: it 
has been proposed that the Sun may have 
uncharacteristically high elemental abun- 
dances (e.g., Mathis 1996), and that the 
lower mean abundances in B stars may 
be more appropriate (Cunha & Lambert 
1992, 1994; Holmgren et al. 1990). The 
recent downward revision of the solar C, 
N, and O abundances may bring the inter- 
stellar standard back closer to solar (Sofia 
k Meyer 2001; Allende Prieto et al. 2001, 
2002). Gas in the Galactic halo in par- 
ticular has less than solar abundances, al- 
though if the hot gas being observed is re- 
cently ejected from the disk, it is less clear 
what adjustment needs to be made. If el- 
ements are depleted onto dust grains, this 
will most significantly affect N(Si IV). It 
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will also somewhat decrease N(C IV), and 
N(C IV)/N(0 VI) (see discussion in Slavin 
& Cox 1992, in the specific context of con- 
ductive interface models). 

3. Velocity-Resolved Models of Li- 
like Ions 

The Li-like ions C IV, N V, and O VI are 
produced in nonequilibrium physical situ- 
ations, such as shocks, conductive inter- 
faces, and rapidly cooling gas (§ 2). The 
different ions are generally produced at dif- 
ferent flow velocities in the structure, so 
a viewer looking through such a structure 
will see trends in the ion column density 
ratios function of line-of-sight veloc- 
ity. The ratio-velocity signature in rele- 
vant physical situations must be consid- 
ered to properly interpret trends in the 
data. Such models are presented in the 
following pages. For conductive interfaces, 
we took existing models and calculated 
the velocity signature that would be ob- 
served along a sightline passing through 
the interface. For radiative shocks, we ran 
the shock code of J. Raymond and calcu- 
lated the velocity signature. For super- 
nova remnants, we consider the velocity- 
resolved results previously presented by 
Shelton (1998). Finally, for a Galactic 
fountain, we integrated the steady-state 
equations published by Houck & Bregman 
(1990), coupled with the recombination 
code of Benjamin et al. (2001), to deter- 
mine the ionization structure as a function 
of line-of-sight velocity. 

3.1. Conductive Interface 

One of the oldest interface models con- 
taining the Li-like ions is that of a con- 
ductive interface. Figure 2 shows the 



log of column density ratios as extracted 
from the models of Bohringcr & Hartquist 
(1987). We created this figure by extract- 
ing the Li-like ion fraction, temperature, 
outflow velocity, and pressure as a function 
of radius from their plots and calculating 
the column density of the Li-like ions at 
each radius. Then, we integrated the col- 
umn density over radius, with a parcel of 
gas at a certain radius modeled in the ab- 
sorption spectrum as a Gaussian velocity 
component centered at the outflow veloc- 
ity at that radius and having the thermal 
width of the gas at that radius. 

Borkowski et al. (1990) present more 
sophisticated conductive interface models, 
with time evolution and a magnetic fleld. 
Although the authors do not provide spa- 
tially resolved ionization information, they 
note that the thickness of an ionized layer 
increases with ionization stage, with G IV 
occupying ~0.3 pc, N V 1-3 pc, and O VI 
~8 pc. During the initial evaporative 
phase, the velocity structure shown in their 
Figure 4 appears similar to the models of 
Bohringer & Hartquist (1987), and would 
thus also show a decrease of ionization 
stage with increasing outward velocity. It 
is unfortunately not possible to verify this 
from the limited information presented in 
Bohringer & Hartquist (1987). During the 
steady-state and condensing phases, there 
is flow back from the hot medium, and thus 
for the longest-lasting phase of the model, 
increasing ionization towards more nega- 
tive velocities would apparently still be ob- 
served. From the point of view of an ob- 
server in the cool cloud, N(G IV)/N(0 VI) 
and N(N V)/N(0 VI) would increase to- 
wards more positive velocities as in the 
Bohringer & Hartquist (1987) models, but 
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Fig. 2. — Column-density ratio as a function 
of outflow velocity calculated from the model of 
Bohringer & Hartquist (1987) of an evaporating 
cold cloud in a hot medium. An observer in the 
cool cloud would sec this ion ratio spectrum, with 
positive velocities moving away from the observer. 
The thickness of the N V/0 VI line scales with 
the density at that point in the flow. Of the sight- 
lines in this study with all three ions, some have 
N V/0 VI profiles that can be explained by a 
conductive interface, but the models predict lower 
C IV/0 VI than observed. 



the absorption would be seen at overall 
negative velocities. Naturally the situa- 
tion is reversed for an observer in the hot 
medium. 

An interesting note is that the models 
of Borkowski et al. (1990) predict increas- 
ing linewidths with ionization potential. 
Linewidths are primarily thermal, and the 
mean temperature at which a certain ion 
is found increases with ionization poten- 
tial. The effect is much more pronouced 
(~25% increase in linewidth) from C IV 
to N V than from N V to O VI (-4% in- 
crease). The observations of Savage ct al. 
(1995) and Sembach & Savage (1992) note 
an increase in line width between C IV and 
N V along sightlines that pass through the 
disk and halo of the Galaxy. New observa- 
tions (paper II) also show a hint of down- 
turn in the N(C IV)/N(0 VI) ratios in the 
line wings, consistent with the trend seem 
by other investigators for narrower C IV 
linewidths than other Li-like ions. No such 
trend is seen in the N(N V)/N(0 VI) ra- 
tio, but in conductive interfaces similar to 
the models of Borkowski et al. (1990) and 
in gas in coUisional ionization equilibrium, 
the expected linewidth difference is much 
smaller between N V and O VI than be- 
tween C IV and O VI. 

3.2. Radiative Shock 

Shock waves are ubiquitous in the in- 
terstellar medium. Of particular interest 
here are the possible shock at the edge of 
the local bubble and shocks in a Galac- 
tic fountain flow. The latter could arise 
on the way up, either in the tops of su- 
pershells that arc bursting out of the disk, 
oblique shocks in more complicated struc- 
tures which are breaking out, or on the 
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way down, with condensed clouds reshock- 
ing against the diffuse halo as they return 
to the disk. There is evidence of O VI in 
high velocity clouds (HVCs), which could 
result from shocking as they fall through 
the halo. However, if the halo is close 
to the virial temperature, then the Mach 
number of the HVCs is only Ms ~1.2-1.5, 
marginally high enough to shock ionize up 
to O VI. 

Figure 3 shows the structure of a 200 km s"^ 
radiative shock, which we calculated us- 
ing the Raymond shock code (J. Ray- 
mond, private communication 1998, Harti- 
gan et al. 1987). Gas immediately behind 
the shock is highly underionized compared 
to its high temperature. Cooling is slow 
above a few million degrees and the gas 
ionizes up. Further behind the shock rapid 
cooling takes place and the gas recombines 
back through the various stages. There 
is thus a secondary rise in Li-like ions in 
the cooling zone, the degree of which de- 
pends on how much of the He-like ions were 
produced, which depends on the shock 
velocity and the metallicity of the gas. 
The details of the ionization- velocity struc- 
ture seen by an observer looking through 
a shock may change depending on these 
parameters, but the general trend is the 
same, because the bulk of the observed Li- 
like ions are are seen as the gas is ionizing 
up, and C IV, N V, and O VI are ionized se- 
quentially as the flow velocity slows down. 
Figure 4 illustrates the ion column density 
ratio as a function of rest-frame flow ve- 
locity, i.e. the hne-of-sight velocity seen 
looking through a shock moving towards 
the observer. As before, the thickness of 
the line scales with the density of N V 
and O VI present along the line of sight. 



so although there are parts of the shock 
in the cooling zone in which the column 
density ratio N(N V)/N(0 VI) increases 
with velocity, in the part of the shock with 
the largest column of those ions, the ra- 
tio is decreasing with velocity. A negative 
slope of N(N V) /N(0 VI) with increasingly 
positive velocity would be observed in a 
sightline dominated by a a shock moving 
towards towards the observer, for example 
on the front faces of Galactic fountain or 
high velocity clouds falling back to the disk 
and being shocked against the diffuse halo 
gas. 

3.3. Supernova Remnant Shell 

Another structure which may be partic- 
ularly relevant to observations of Li-like 
ions in the Galactic halo is a supernova 
remnant shell. The tops of superbubble 
shells may be launched into the halo as 
partially coherent clouds, and if we are 
living inside a superbubble, there may be 
remnants of the top of our shell above us, 
through which many of the sightlines in 
this study would pass. 

Shelton (1998) modeled evolving super- 
nova remnant (SNR) bubbles expanding 
in ambient conditions typical of the lower 
halo. In particular, she presented the 
velocity-resolved column densities of Li- 
like ions (her Figure 11). Those values and 
the column density ratios are shown in Fig- 
ure 5. Clearly, the ion-ratio proflle depends 
on the stage of evolution of the SNR. At 
very early stages, the interior of the rem- 
nant has not yet ionized past the Li-like 
stages, so these species are found both be- 
hind the shock and in the interior, lead- 
ing to the double-peaked proflle in column 
density. Later, the ions are found in the 
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distance (lO'^cm 



Fig. 3. — Structure of a 200 km s"^ interstellar 
shock, calculated using the Raymond shock code. 
The first plot shows the Li- like ionization fractions 
of C, N, and O, the temperature, and the flow ve- 
locity as a function of distance behind the shock. 
The second plot shows the ionization fractions of 
all stages for oxygen. Intermediate and high ion- 
ization stages (e.g., O V here) are present imme- 
diately behind the shock as the gas ionizes up, and 
farther back in the recombination zone, as the gas 
cools back through ~10^ K. 
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Fig. 4. — Ion column density ratios as a function 

of flow velocity in a 200 km s'^ radiative shock. 
The thickness of the line scales with the density of 
N V and O VI present along the line of sight. The 
first plot shows the column density ratio which 
would be measured by an observer capable of mea- 
suring the volume density in situ. The second plot 
shows what would be measured by an observer 
using absorption spectroscopy through the shock. 
The gas at each point in the flow is convolved with 
the thermal width at that point in the flow, greatly 
smoothing the observed structure. 
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cooling but overionized gas at the periph- 
ery of the remnant. Finally, as the aged 
remnant cools, O VI and N V are found 
mostly in the interior, and C IV near the 
outskirts. 

The profiles of aged remnants are most 
relevant to the possibility of our looking 
out through our own superbubble shell, 
and these profiles show a decrease in 
N(N V)/N(0 VI) with increasing velocity. 
Additionally, if there is a rising shell above 
our superbubble, it is probably fragmented 
owing to Rayleigh- Taylor instabilities as 
the superbubble blowout expands through 
the density gradient out of the Galactic 
disk. Fragmentation and further distur- 
bance by Kelvin-Helmholtz instability, as 
the shell passes through more stationary 
diffuse gas, would confuse any velocity- 
ionization signature. The nature of our 
local bubble would be better tested with 
spectroscopy towards stars just outside the 
local bubble, in which the gas at the bubble 
interface could be observed independently 
from other gas which may be involved in 
different physical processes. The theoret- 
ical ion-ratio profiles of remnants at dif- 
ferent ages could be very useful for such 
local bubble observations, or other obser- 
vations in the disk of the Galaxy which 
pass through known remnants. If the ge- 
ometry is well-enough understood, obser- 
vations of the column density ratio could 
determine the age of the remnant. 

3.4. Galactic Fountain 

The nature of a Galactic fountain flow 
is determined by how the quickly gas cools 
and recombines as it flows upwards out of 
the disk. Quantitatively, the comparison is 
between the cooling time Tc and the sound 



crossing time Tg, which is roughly the time 
to establish hydrostatic equilibrium: 



and 

nkT 

" (7-lKni^A(r)' 

where 7 is the adiabatic exponent, /i is 
the atomic weight, and A(T) is the cool- 
ing function (e.g., Houck & Bregman 1990, 
who discuss the nature of the fountain 
flow for the different values of Ts/tc). If 
Tc S> Ts then a static system can be estab- 
lished, but it will be unstable to convective 
overturn and thermal instability (as noted 
by Bregman 1980, the cooling rate of the 
static structure is highest at the top of the 
halo, so clouds will condense out from the 
top down and be formed with small neg- 
ative velocities). If ^ then a su- 
personic ffow will occur. As the gas cools 
rapidly it is thermally unstable and will 
condense into clouds which fall out of the 
ffow. The scale height of the gas is then set 
by the cooling time H ~ r^v ~ r^Cs, where 
Cg is the sound speed at the base of the 
fountain, and clouds can be formed with 
positive (upward) velocities. For Tc ~ r^, 
transsonic fountain ffows exist with clouds 
forming near the adiabatic density scale 
height. 

Modeling of the velocity-ionization struc- 
ture of a Galactic fountain will focus here 
on the transsonic solution (r^ ~ Tc) of 
Houck & Bregman (1990). (As discussed 
below, ffow solutions in which the cooling 
and equilibrium times are very different 
are not as well-suited to analytic treat- 
ment.) Figure 6 shows the structure of the 
Houck & Bregman (1990) fountain flow, an 
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integration of the one-dimensional steady- 
state fluid equations, which pass through a 
sonic point. The equations have the same 
topology as a standard solar/stellar wind 
(Parker 1958). The solution stagnates at 
the top, forming an unphysical dense layer, 
but in this zone of rapid cooling clouds will 
condense out with near-zero or small posi- 
tive vertical velocities and fall back to the 
disk. 

Figure 7 shows the oxygen ioniza- 
tion fractions in a nonequilibrium cooling 
fountain flow, which we calculated using 
the nonequilibrium recombination code of 
Benjamin et al. (2001). We integrated the 
recombination rate through the tempera- 
ture and electron density given as a func- 
tion of time and vertical velocity by the 
Houck & Brcgman (1990) transsonic foun- 
tain. It would be slightly more consistent 
to integrate the fluid equations with Ben- 
jamin et al. (2001) nonequilibrium cooling 
instead of the Edgar (1986) nonequilibrium 
coohng used by Houck & Bregman (1990), 
but the differences in the cooling curves 
are small enough for the combination to 
be accurate (R. Benjamin, private com- 
munication 2001). O VII cannot recom- 
bine fast enough as the flow cools rapidly 
through its coUisional ionization equilib- 
rium temperature,~(5-10)xl0^ K. Signif- 
icant O VII remains to low temperatures 
(<10^ K if the density is low enough), until 
dielectronic recombination allows a recom- 
bination cascade. The result is that not 
only does O VI persist to low temperatures 
because it is "frozen in" (it cannot recom- 
bine fast enough to keep up with cooling), 
but there can be a rise in the O VI frac- 
tion as the He-like O VII recombines. This 
happens for the other Li-like ions N V and 



^ ^ Early epoch (10* yr) ShjII formation (2. 5x 1 O' yr) K 
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Fig. 5. — Column densities and ratios as a func- 
tion of velocity for a sightlinc through the center of 
a supernova remnant, at various stages of the bub- 
ble evolution (Shelton (1998)). The profiles are 
symmetric about zero velocity, and only positive 
velocities arc shown. Note that the velocity scales 
change significantly as the shock decelerates. 




Fig. 6. — Houck & Bregman (1990) transsonic 
fountain. The vertical dotted line marks the sonic 
point. 
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C IV for the same reasons. The magnitude 
of this effect is dependent on dielectronic 
recombination rates which have undergone 
continual revision with time (Benjamin 
et aL 2001; Shapiro & Moore 1976; Kafatos 
1973), and makes the modehng of Li-hke 
ions in cooled (T^IC^ K) clouds uncertain 
(see discussion and a crude model below). 

Figure 8 shows the ion column density 
ratio calculated by coupling the Benjamin 
et al. (2001) reionization with the Houck & 
Bregman (1990) transsonic fountain. The 
parameters for this particular model are 
set by choosing the fountain base temper- 
ature to be 6x10^ K. The two branches 
in the first plot are a result of the veloc- 
ity of the fountain accelerating, then de- 
celerating again after passing through the 
sonic point. The horizontal upper branch 
is from the rapidly cooled gas at a few 
10^ K, which will condense into clouds and 
fall out of the flow. The column density 
of Li-like ions at this point in the flow is 
highly uncertain, and the column plotted 
here is probably an overestimate. (R. Ben- 
jamin, private communication 2001). The 
width of the line in the plot scales with 
the column density of N V and O VI at 
that point in the flow. Thus, an observa- 
tion through the fountain is dominated by 
the rising, cooling gas. The second plot 
shows the ion column density ratio which 
would be observed using absorption spec- 
troscopy. Each parcel of gas in the flow 
was convolved with the thermal width at 
that point in the flow to create an ob- 
served absorption profile, and the ion ra- 
tio is calculated from that observed pro- 
file. Since the transsonic fountain acceler- 
ates, N(N V)/N(0 VI) increases with in- 
creasingly positive velocity. However, the 
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time (1 0^ years) 

Fig. 7. — Oxygen ionization fractions in a 

noncquilibrium cooling fountain flow, calcnlatcd 
using the analytical approximation to the hydro- 
dynamic equations of a transsonic fountain of 
Houck & Bregman (1990) coupled to the ioniza- 
tion code of Benjamin et al. (2001). O VI persists 
to low temperatures because the gas cools much 
faster than it can rocombinc. The O VI ioniza- 
tion fraction can even rise at low temperatures due 
to the much-delayed recombination of O VII (see 
text). 
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trend is very weak due to the small velocity 
gradient in this fountain model, and it is 
not clear that the transsonic smooth foun- 
tain is a realistic description of the patchy 
Galactic halo. 

It has been noted (de Avillez 2000) that 
there have been two types of Galactic foun- 
tain models in the literature, those with 
somewhat smoothly distributed rising gas 
over the disk, as used here, and models 
that focus on the evolution of superbubbles 
as the fountain drivers, irregularly spaced 
in time and galactic location (e.g., Nor- 
man & Ikeuchi 1989). Although the lat- 
ter model is driven differently and rises 
to larger Galactic altitude than that of 
Houck & Bregman (1990), the cooling and 
equilibrium times are still close, the ex- 
pected spread of velocity in the rising gas 
is modest, and so the observed ion col- 
umn density ratio as a function of veloc- 
ity is probably not qualitatively different 
from the model considered here. What 
is needed is an improvement over either 
type of simple model, i.e. full multidi- 
mensional time-dependent models of the 
fountain with nonequilibrium cooling and 
nonequilibrium reionization. 

Galactic fountain flows in which the 
cooling and equilibrium times are very dif- 
ferent are less well-suited to an analytic 
treatment because, as discussed above, it 
is difficult to accurately follow the recom- 
bination state of gas which is highly ovcri- 
onizcd but which has cooled to ~10'^ K and 
thermally condensed. In a supersonic foun- 
tain (tc < Tg) clouds condense quickly and 
observations are dominated by this uncer- 
tain cool phase. In a highly subsonic foun- 
tain (tc ^ Tg) there is no velocity signa- 
ture because the gas rises and establishes 
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Fig. 8. — Ion column density ratios as a func- 
tion of velocity for a Houck & Bregman (1990) 
transsonic fountain coupled to Benjamin et al. 
(2001) nonequilibrium reionization. The thickness 
of the line scales with the N V and O VI den- 
sity along the line of sight, so the observed ab- 
sorption profile is weighted towards the ratio with 
the thickest line. For clarity, the column density 
weighting is only shown for log[N(N V) /N(0 VI)] . 
The first plot shows the ion column density along 
the flow as it would be measured by an observer 
in situ capable of measuring the ion volume den- 
sity at a point in the flow. The second plot shows 
what the observed ion column density ratio would 
be for an observed using absorption spectroscopy 
as was done in this study. The ion column den- 
sity at each point in the flow is convolved with the 
thermal width of the gas at the temperature at 
that point in the flow to create an absorption pro- 
file, and the trend is thus significantly smoothed 
out. 
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a quasistatic configuration before cooling 
through the Li-hke ions. Any difference in 
ionization function of velocity again 
occurs in the cold clouds. 

Nevertheless, to estimate the behavior 
of recombination in cool clouds, Figure 9 
shows the ion column density ratio as a 
function of time for a crude model. We 
consider a cloud formed by thermal insta- 
bility in a smooth fountain flow at 10"^ K 
with N V, N VI, O VI, and O VII ioniza- 
tion fractions of 0.06, 0.01, 0.02, and 0.005 
respectively. These values are fairly un- 
certain, and were taken from Shapiro & 
Moore (1976) with consideration of Ben- 
jamin et al. (2001) and Sutherland & Do- 
pita (1993) and the known limitations of 
these models (R. Benjamin, private com- 
munication 2001). The cloud cools slowly 
at a few xlO^ K, remaining in pressure 
equilibrium with the rising smooth hot 
halo, and thus the density increases, but 
not by a huge factor. At the same time, 
hydrogen recombines, so the electron den- 
sity must decrease. Considering the other 
large uncertainties of this model, the sim- 
ple approximation is adequate to allow hy- 
drogen to recombine down to the hydro- 
gen ionization fractions Xh ~ 0.1 — 0.5 
observed in diffuse halo clouds (Spitzer 
& Fitzpatrick 1993) and modeled in pho- 
toionization equilibrium with ionizing ra- 
diation escaping from the Galactic disk 
(Bland-Hawthorn & Maloney 1999) or old 
supernova remnants (Slavin et al. 2000). 
Total (radiative and dielectronic) recombi- 
nation rates were mostly taken from (Na- 
har et al. 2000; Nahar 1999; Nahar & Prad- 
han 1997; Schippers et al. 2001), which 
shows significant differences from previous 
work. 
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Fig. 9. — Ion column density ratio 

log[N(N V)/N(0 VI)] as a function of time 
in recombining cool gas. Initial ionization 
fractions are indicative of what might be found 
in gas that has rapidly cooled to a few 10'* K, 
is highly overionized, and has condensed into a 
cloud by thermal instability. The exact behavior 
(for example, the magnitude of the rise relative 
to the drop in ion ratio) depends on the choice 
of uncertain recombination coefficients and even 
more uncertain initial ionization states. 
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The ion ratio N(N V)/N(0 VI) first de- 
creases with time, as the N V total recom- 
bination rate is larger than that for O VI. 
This fact was not true for the older rates 
of ShuU & Van Steenberg (1982). The col- 
umn density ratio then rises as additional 
N V and O VI rccombine down from N VI 
and O VII (the recombination rate for the 
latter is larger). As elsewhere, the thick- 
ness of the line in Figure 9 scales with 
the local density of N V and O VI. The 
behavior of the ion ratio with time (how 
much it rises at late times, etc.) varies 
with the choice of initial ionization frac- 
tions, but the recombination time in the 
cool cloud is always short compared to the 
return time for the ballistically traveling 
cloud: the cloud velocity in the Galactic 
gravitational field {g ^10"*^ cm s"^) only 
changes by 3 km s"^ in the 10^ year re- 
combination time. Thus any change in the 
ionization ratio function of velocity 
determined by the distribution of vertical 
velocities of clouds at their formation time, 
not by a slow change in ionization state as 
the clouds decelerate. The height and ver- 
tical velocity of cloud formation depends 
on the initial temperature of the fountain 
as discussed above, and simulations show 
a fairly large spread for a given fountain 
flow (Bregman 1980). Although the gen- 
eral assertion can be made that clouds re- 
combine as they travel the fountain flow, 
no conclusions can be drawn about an ob- 
servable velocity-ionization signature from 
this type of simple considerations of single 
clouds. What is needed are full nonequilib- 
rium 3-d simulations of fountain dynamics. 



3.5. Other Halo Models 

Although the most consistent model 
of the halo is of a dynamic flow of gas 
heated and ionized by (thermal) collisions, 
it is worth considering some of the other 
physics that could be important and how 
that might relate to observed Li-like ion 
column density ratios. 

Magnetic fields play a large, but usu- 
ally neglected, role in interstellar dynam- 
ics. For example, the magnetic pressure in 
the Galactic disk is about the same as the 
thermal pressure. Synchrotron emission 
from our Galactic halo and the halos of 
other Galaxies show that the scale height 
of the magnetic field is large (Han & Qiao 
1994). Dissipation of magnetic energy is a 
viable mechanism for heating gas to mil- 
lions of degrees in the disk (Tanuma et al. 
2001; Hartquist & MorfiU 1984), moder- 
ately ionizing and heating diffuse gas in 
halo (Birk et al. 1998), and even in produc- 
ing the X-ray emission of boundary regions 
of high velocity clouds passing through 
the halo (Zimmer et al. 1997). Raymond 
(1992) presents a quasistatic halo model in 
which the gas is heated by magnetic re- 
connection. The model can produce rea- 
sonable column densities of high and mod- 
erately ionized species, and reasonable in- 
tensities for emission lines of those species. 
It may have some difficulty fitting both 
the highly ionized species and the less ion- 
ized species, because the lower energy re- 
connection mechanism which is invoked re- 
quires possibly more magnetic power than 
is available in the Galaxy. In addition, 
this model would have no kinematic signa- 
ture, being quasistatic, and the stability of 
such a model is not discussed in Raymond 
(1992). 
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Cosmic rays are another sometimes ne- 
glected but important part of the interstel- 
lar medium. Again, the cosmic ray pres- 
sure in the plane is comparable to the ther- 
mal and magnetic pressures. Boulares & 
Cox (1990) presented a simple static halo 
supported by cosmic rays (CR), and ar- 
gued that the stability times arc of the or- 
der of 10^ years. More recent work calls 
into question the stability of CR supported 
haloes, especially in the presence of CR dif- 
fusion which makes CR supported struc- 
tures more Rayleigh- Taylor unstable (Bre- 
itschwerdt et al. 1993). A static CR sup- 
ported halo would also not have variation 
of ion column density ratios with line-of- 
sight velocity. Perhaps more likely is a CR- 
driven wind, which may operate in paral- 
lel with a Galactic fountain, driving ma- 
terial out in regions of open magnetic field 
topology, while lofted material in regions of 
closed magnetic field falls can fall back to 
the disk (Breitschwerdt et al. 1991). These 
authors concentrate on the CR-driven out- 
flow at hundreds of kiloparsecs from the 
disk, but their model does have monoton- 
ically increasing velocity as gas rises and 
cools, entrained in the CR wind. The gra- 
dient of ~4 km s"^ kpc"^ is almost identical 
to that in the lower part of the Houck & 
Bregman (1990) transsonic fountain, and 
thus would have a very similar trend of ion 
column density ratio with velocity. There 
is the possibility that damping of hydro- 
magnetic waves in the outflow can trans- 
fer energy to the gas; the recombination 
of the gas could be thus delayed to higher 
altitudes, where the velocity gradient is 
smaller, and perhaps a steeper change in 
ion column density ratio as a function of 
velocity could be produced. To fully un- 
derstand the Li-like ions in this CR wind. 



more specific discussion of ionization is re- 
quired than is provided in Breitschwerdt 
et al. (1991). 

Finally, one can consider that the halo 
gas is partially photoionized. It has been 
argued that photoionization is not negligi- 
ble in the halo (Slavin et al. 2000; Bland- 
Hawthorn & Maloney 1999), but it is diffi- 
cult to explain species with large ionization 
potentials by photoionization, namely N V 
and O VI. This can be easily shown by cal- 
culating the path lengths required to pro- 
duce the observed N(N V)/N(0 VI) col- 
umn density ratio and simultaneously the 
total column, or the N(C IV)/N(0 VI) ra- 
tio. The conclusion is similar for the entire 
data set: ion column density ratios require 
high ionization parameters and low space 
densities, but substantial column densities 
then require path lengths of tens or even 
hundreds of kiloparsecs, which is unreason- 
able to fit in the Galactic halo. In addition, 
the two ion ratios N(C IV)/N(0 VI) and 
N(N V)/N(0 VI) cannot simuhaneously 
fit a photoionization model in most cases. 
The ionization parameter U — n^/nu {n^ 
is the number density of hydrogen-ionizing 
photons) derived from N(C IV)/N(0 VI) is 
about 0.2 dex lower than that derived from 
N(N V)/N(0 VI), indicative of a different 
(coUisional) ionization source. 

4. Conclusion and Summary 

We have considered the diagnostic power 
of velocity-resolved column density ratios 
in understanding the Galactic halo. Col- 
umn density ratios of Li-like ions in the 
Galaxy are useful to diagnose the physi- 
cal formation mechanism of the gas and 
to study the interstellar gas cycle, and 
a survey of these ions can reveal general 
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trends. In Paper II, we present a survey of 
sightlines observed with FUSE and HST, 
in which the distribution of N(N V) and 
N(0 VI) in the halo does not appear to fa- 
vor a dominant physical production mech- 
anism. 

Here, we have presented models of inter- 
faces and cooling nonequilibrium gas, fo- 
cusing on the velocity-resolved N(N V) /N(0 VI) 
signatures. In particular, we consider con- 
ductive interfaces, radiative shocks and 
supernova remnants, and cooling gas in 
a Galactic fountain. Typical ion column 
density ratios are summarized in Table 3, 
along with a typical observed slope ob- 
served in FUSE and HST data. One im- 
portant type of hot/cool gas interface for 
which the ionization- velocity signature has 
not been discussed is the turbulent mixing 
layer (Slavin et al. 1993). The structure 
of these Li-like ion producing structures 
is uncertain. It is likely that there is a 
turbulent cascade of eddies which would 
wash out any velocity-ionization signature, 
but there may possibly be large Kelvin- 
Helmholtz rolls with some coherent veloc- 
ity signature. 

The observable velocity-ionization trends 
are weak, because even very strong trends 
are washed out by the large thermal width 
of the gas at different parts of the flow. 
These trends could be further comphcated 
by thermal instabilities, which likely oc- 
cur in cooling gas such as in a fountain 
flow, and can also affect shocks and super- 
nova shells. Fragementation of the gas into 
parcels with complex density and velocity 
structure, as seen in the time-dependent 
shock model of Sutherland et al. (2003) 
could complicate the signatures described 
here. 



Additional confusion can result when 
long sightlines pass through multiple struc- 
tures. As we show in Paper 11, the dis- 
persion of N(N V)/N(0 VI), both inte- 
grated and velocity-resolved, indicates that 
no single production scenario known to 
date can completely explain the Galac- 
tic halo. To truly understand the physi- 
cal production of Li-like ions in the halo, 
one needs to analyze gas in localized ar- 
eas of physical space, rather than velocity 
space. Absorption spectroscopy towards 
many halo stars with close angular separa- 
tion and different distances could help to 
isolate gas at a specific altitude. Similarly, 
the gas above known superbubble shells or 
chimneys could be isolated. These observa- 
tions have a greater chance of distinguish- 
ing between models of hot gas production 
than observations along long lines of sight. 
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Table 3 
Summary of models 



Model 


Slope of Ratio log[N(N V)/N(0 VI)] 




Radiative shock 


~-0.0015 dex (km s" 


::>:: - 


Mature SNR shell (Shehon 1998) 


>-0.02 dex (kin s' 




Ilouck & Br(-guiaii (1990) iouulaiii 


~+ 0.002 (lex (km s" 




Indicative observed slope 


-0.0032+/-0.0022(r)+/-0.0014(s) dex (km s" 





Note. — Models 
through the structu 
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